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Abstract 


A realistic nuclear mean-field hamiltonian with pairing has been di¬ 
agonalized using Fock space representation that allows for nearly exact 
treatment of the problem. Calculations were performed for all the even- 
even nuclei with Z G (20,100), whose pairing gaps were possible to extract 
out of the experimental masses. The optimal values of the pairing strength 
constants for the protons and neutrons have been found. 

PACS: 21.30.Fe, 21.60.-n, 71.10.Li 


In the large scale microscopic calculations of the nuclear total-energy surfaces 
the method of Strutinsky plays an important role allowing for fast and fully au¬ 
tomatic calculations related to the equilibrium deformations, shape coexistence, 
fission probabilities and many other mechanisms and phenomena. In the related 
formalism the calculation of shell- and pairing-energies plays a decisive role, the 
latter obtained so far with the help of the Bogolyubov transformation and the 
associated Bardeen-Cooper-Schrieffer approximate method. The new method 
proposed in Ref. [|lj allows to obtain the exact (in some cases nearly exact) so¬ 
lutions of the pairing problem using realistic hamiltonians - in particular those 
with the state-dependent pairing hamiltonian. The new method is based on 
the direct solution of the many-body problem in Fock space; it employs tech¬ 
niques similar to those used in the nuclear shell-model, including the Lanchos 
diagonalisation scheme. The hamiltonian in question has the form: 
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where H m f denotes any mean-field hamiltonian e.g. the one with deformed 
Woods-Saxon potential and matrix G a j 3 , in general non-diagonal in its in¬ 
dices, defines the state-dependent pairing-hamiltonian. Owing to three exact 
symmetries obeyed by hamiltonians of the above general form, cf. Ref. pj for 
details, the corresponding matrix written down using a Fock space basis can be 
block-diagonalized analytically, thus reducing the problem of diagonalisation of 
huge-size matrices to much smaller ones that can be treated easily with the help 
of the Lanchos methods. 

In this article we report on the results of an introductory large-scale test 
that consists in fitting the coupling constants of the simplest (monopole-pairing) 
version of pairing hamiltonian in (Q), i.e. G «,/3 = G5 a ^. Pairing calculations 
were performed within Fock space defined by 24 particles placed on 24 double¬ 
degenerate single-particle energy levels in the ’pairing window’. The corre¬ 
sponding full hamiltonian-matrix has the size of Nham = 32 247 603 683100, 
while the sizes of hamiltonian-blocks, after applying the formalism of Ref. [jl|, 
are: N s —q = 2 704 156 in the seniority-zero and N s —2 = 705 432 in seniority-two 
blocks. In both cases we have applied a basis cut-off reducing the sizes of the 
effectively diagonalized matrices to fV s ' =0 = 27 703 and fV s ' =2 = 26 263. We 
have verified by comparison with the results of the exactly soluble method of 
Richardson that the basis cut-off that reduces the matrix sizes by roughly two 
orders of magnitude induces errors of ~ 2% only. 

We have diagonalized the Fock-space hamiltonian-matrix corresponding to 
( 0 ) and we obtained the ground-state energy (seniority-zero) and the lowest- 
energy seniority-two states. The difference between the two energies has been 
interpreted as corresponding approximately to twice the ’pairing gap’, the latter 
obtained from the mass-difference expression 


A^(N) = -y(B(N) + B(N + 2) - 2 B(N + 1)), 


( 2 ) 


as discussed recently in detail in Ref. I- Above, 7 t/v = (—1)^, B(N) are (neg¬ 
ative) binding energies of nuclei, for the fixed ^-number. To obtain the proton 
’pairing gap’ one has to replace N with Z and fix the neutron number N. The 
resulting proton and neutron pairing gaps for even-even nuclei are illustrated 
in Fig. jl], top, in function of the mass number A; the A ^-values with the 
experimental errors exceeding 250 keV were not taken into account. 

The best average fit in terms of the a/\/~A one-parameter dependence for 
all the nuclei is = 10.85/\/A and Ap 3 ^ = 10.02 /y/~A (solid lines in Fig.llj, 
top). The fitted a-values are smaller than those obtained in [3] from the liquid 
drop model formula, the corresponding curves denoted in Fig.l with the dashed 
lines. One can also use alternative one-parameter expressions similar to those 
in Ref. 101; the results of the fit are 


or 


A n = 23 /{VAf and A p = 21 /(VA) s , (3) 

A n = 4.66/A 1 / 3 and A p = 4.18/A 1 / 3 , (4) 

but the fit precision remains similar to that with the a/y/A-type dependence. 
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Figure 1: Experimental A p , protons, frame A), and A„ , neutrons, frame 
B), together with the average-fit curves of the form a/\/~A (top). The usually 
quoted average behavior \2/y/~A is also illustrated. Frames C) to F) represent 
the results of the pairing constant-G fitting procedure (see the text and Table 
|T|) . Experimental masses are from Ref. ji|. 


To avoid the undesired type of variation in pairing delta, characterististic for 
spherical (especially doubly magic) nuclei, the nuclei with available experimental 
masses have been arbitrarily divided into 5 regions: (I) with Z G (32, 38) and 
N G (32,44); (II) with Z G (40,46) and N G (56,72); (III) with Z G (54,66) 
and N G (58,76); (IV) with Z G (60,82) and N G (88,104) and (V) with 
Z G (90,100) and N G (142,156); there the doubly-magic nuclei have been 
eliminated. Within those regions the A-values were extracted and, as the next 
step, the G-constants found that reproduce the extracted A values exactly for 
each of the studied nuclei at the calculated in advance equilibrium deformations. 
The mean field hamiltonian used is the same as in Ref. p (see also references 
there). The irregular behavior of A^ 3 ) in function of A suggests that the resulting 
G values will also vary in a relatively irregular fashion and as a consequence we 
have tried ’a few parameter’ fits in terms of A and I = (N — Z)/(N + Z) and 
possibly some powers of those variables: 

G p(u) = ^jriPo +Pil + P 2 I 2 ) where s = 1,1/2, 2/3. (5) 
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In Fig. |T|, frames C) and D), the products G p ( n ) • A are shown in function of the 
isospin factor I while the corresponding average behavior is shown with the help 
of the solid lines. In frames E) and F) the quantities G p ( n ) ■ A 2 / 3 are illustrated. 
The dependence in terms of / obtained here is nearly constant. 


Table 1: Pairing strength constants in terms of the approximating expressions 
for neutrons and protons for 5 regions considered. 


Gp(n) 

a/A 4- bl 

(c + d I)/ A' 2/ J 

(c 1 + d 1 1 + el 2 )/A 213 

Region 

a 

b 

c 

d 

d 

d' 

e 

alln 

24.872 

0.147 

6.517 

-5.921 

5.657 

10.96 

-64.18 

I 

24.230 

-0.001 

5.887 

-1.438 

5.657 

7.43 

-65.63 

II 

27.177 

-0.165 

6.230 

-6.453 

5.657 

0.67 

-21.39 

III 

33.129 

-0.214 

7.097 

-10.298 

5.657 

19.60 

-152.21 

IV 

56.170 

-0.849 

11.457 

-33.690 

5.657 

36.46 

-206.93 

V 

38.548 

-0.144 

6.315 

-6.408 

5.657 

-0.29 

-13.35 

allp 

26.861 

0.331 

7.100 

-2.905 

6.529 

2.65 

-5.51 

I 

24.548 

0.018 

5.934 

4.493 

6.529 

-9.13 

103.58 

II 

45.893 

-0.750 

7.249 

-2.892 

6.529 

1.74 

-3.91 

III 

39.296 

-0.481 

7.178 

-3.700 

6.529 

4.55 

-7.60 

IV 

54.546 

-0.608 

7.289 

-3.327 

6.529 

1.88 

-4.47 

V 

78.763 

-0.635 

8.333 

-4.289 

6.529 

6.18 

-8.67 


Results in Table [l] give the overall r.m.s. deviations that are rather small. 
The variations obtained, based on the experimental data concerning a broad 
range of nuclei, show more structure than the simple parametrisations tested 
can take care of. In particular, strong variations in parameter e from one nuclear 
range to another, cf. column 8 in the Table, indicates that the I 2 fluctuations 
are too rapid to allow deducing any systematic trends in this context. 

The parametrisations summarized in the Table are ’ready to use’ in conjunc¬ 
tion with the Fock-space diagonalisation method of Ref. jlj. 
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